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ABSTRACT 

We present results for the two-point angular correlation function of galaxies, uj(6), to a limiting 
magnitude of r = 26. Our catalogue is constructed from deep imaging using the COSMIC imaging 
spectrograph on the Hale 5-m. The final sample is 97% complete to r = 26.0 yielding ~5700 galaxies 
over a 90.1 sq. arcmin field. Our analysis shows uj(6) for faint galaxies can be parameterised by 
a power law of the form A^O -0 ' 8 , in agreement with the angular clustering statistics of shallower 
catalogues. The derived amplitude, A^, for our catalogue is small, but non-zero. We combine 
this measurement with the latest statistical constraints on faint galaxy redshifts from gravitational 
lensing studies, which imply that the bulk of the r < 26 field galaxies should be at redshifts 
z ~ 1. We show that our derived A u is significantly lower than that predicted from the local 
bright, optically-selected galaxy correlation function using the lensing-determined galaxy redshift 
distribution and modest growth of clustering. However, this simplistic model does not include the 
variation in observed morphological mix as a function of redshift and apparent magnitude in our 
sample. At our faintest limits we reach sufficiently high redshifts that differential X-corrections 
will result in the observed galaxy mix being dominated by star bursting dwarf and low surface 
brightness irregulars, rather than the early- type systems used to define the local bright galaxy 
correlation function. Adopting the correlation function measured locally for these low surface 
brightness galaxies and assuming modest clustering evolution, we obtain reasonable agreement 
between our model and observations. This model, therefore, supports the scenario in which the 
high number density of faint galaxies is produced by normally clustered star forming dwarf galaxies 
at modest redshifts. 

Key words: galaxies: clustering - galaxies: evolution - cosmology: observations - large-scale 
structure of the Universe 



1 INTRODUCTION 

The observed surface density of faint galaxies far exceeds the number predicted from a simple 
extrapolation of the local universe (the no evolution model; NE), irrespective of the cosmological 
geometry. However, to the faintest limits measurable the distribution of galaxies in redshift is a 
close match to that predicted by the NE model. Since at fainter limits galaxies become progres- 
sively bluer, the 'excess' population is primarily caused by an increase in the space density of blue 
low-luminosity galaxies at modest redshifts. The exact nature of this blue excess population is 
currently one of the major questions in observational cosmology and has significant impact on our 
understanding of galaxy formation and evolution. 

Observationally, one of the cheapest statistics to obtain for samples of faint galaxies is the 
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strength of their clustering as seen in projection on the plane of the sky. This is commonly measured 
using the two-point angular correlation function, lo(6), defined by 

dP = n 2 [l + Lu(9)}dtt 1 dQ 2 , (1) 

where n is the mean number density of galaxies and dP is the probability in excess of Poisson of 
observing a pair of galaxies in solid angle elements dQ\ and dQ,2 that are separated by an angle 9 
on the sky. To a given depth, u>(9) is dependent upon the redshift distribution of the sources, N(z), 
and their 3-dimensional real space clustering strength, £(r, z). Thus oj{9) is a sensitive test of both 
the luminosity and clustering evolution of faint galaxies (Koo & Szalay 1984). This fact, combined 
with the advent of efficient large format CCD detectors, has led to a recent flurry of interest in 
using the correlation strength of faint galaxies to study their nature (Neuschaefer, Windhorst & 
Dressier 1991; Efstathiou et al. 1991; Couch, Jurcevic & Boyle 1993; Roche et al. 1993). 

At the bright magnitudes accessible with plate material from Schmidt telescopes (B < 20), u>(0) 
has been shown to follow a power law form: uj{9) = AJ^ with 5 w 0.8 (Groth & Peebles 1977; 
Shanks et al. 1980; Maddox et al. 1990). The value of depends upon the depth of the survey, 
but the derived values are in reasonable agreement when the results from different catalogues are 
rescaled to similar depths. Using plates from 4-m class telescopes it is possible to go somewhat 
fainter (B < 24), while still covering a large area. The studies at this depth are in less agreement, 
but do seem to support a continued 5 pa 0.8 exponent for u>{9) to at least B ~ 23 (Koo & Szalay 
1984; Stevenson et al. 1985 (SSFM)). 

To reach still fainter limits CCDs have to be employed, with a resulting large reduction in 
sky coverage. As mentioned, several groups have attempted to measure oj{9) using CCD imaging 
and we briefly summarise these studies. Neuschaefer, Windhorst & Dressier (NWD, 1991) used g 
imaging over a 720 sq. arcmin area complete to g ~ 24.5 (roughly R ~ 23.5) and claim little change 
in the index, 5, of oj{9) from the value at brighter magnitudes. However, they do see a strong 
decline in the amplitude, A^, to g ~ 24.5. Efstathiou et al. (EBKTG, 1991) used multi-colour data 
on a number of small fields (each about 10 sq. arcmin) to measure u>{9) to B < 26 or R < 25. 
Their sample shows very weak clustering at their faintest limits and from this they conclude that 
a large fraction of the faint galaxies in their sample belong to a weakly clustered population which 
is not seen today. By far the largest study is that of Couch, Jurcevic & Boyle (CSB, 1993) who 
analysed a 4 sq. degree area to R < 23 and, as NWD, they found no change in 5 with depth, while 
A^ declined strongly. They conclude that their observations are inconsistent with linear growth in 
a standard CDM cosmogony. Finally, Roche et al. (RSMF, 1993) analysed a total of roughly 300 
sq. arcmin in B and R, with the R sample limited at R < 23.5. Adopting a 9~°- s power law for 
oo(9), they studied the variation of A^ with depth and their analysis shows the, by now familiar, 
decline in A u at fainter magnitudes. They interpret their result as support for a model in which 
the faint galaxies have a clustering amplitude similar to local bright galaxies, but are distributed 
across a wide redshift range: 1 z <^ 3. 

The recent study of Bernstein et al. (1994) used both plate material and statistical information 
on the redshift distribution of galaxies to B ~ 22 to estimate £(r) on small scales at two different 
epochs: z ~ 0.18 and z ~ 0.27. Their results are significantly lower than predictions of models 
extrapolated from local bright galaxy clustering. However, their observations are consistent with the 
assumptions of modest clustering growth and that the bluest 60-70% of their sample have present- 
day clustering similar to IRAS-selected samples, which are intrinsically more weakly clustered 
than bright optically-selected samples (Saunders, Rowan-Robinson & Lawrence 1992, SRRL; Fisher 
et al. 1994). 

A more direct approach to study the evolution of the two-point spatial correlation function, £(r), 
uses the extensive field redshift samples which are gradually becoming available. From the large 
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Broadhurst, Ellis & Colless AAT redshift survey, Cole et al. (1994) analysed both the evolution in 
£(r) and possible differences between the clustering of strongly and weakly star-forming galaxies 
(as indicated by their [On] emission). They detected no evolution in £(r) out to z ~ 0.6 and no 
significant difference between the clustering strengths of the [On]-weak and [On]-strong populations. 
However, at the limit of this study {B ~ 22) neither the apparent surface density of galaxies nor 
their clustering amplitude depart significantly from the NE prediction and, thus, the measurements 
are not particularly sensitive to any excess population. 

The CCD studies of the angular correlation function of faint galaxies, while in broad agreement 
about the strong decline of the clustering strength of galaxies at fainter apparent magnitudes, differ 
in the interpretation of this result. Their inconclusive nature has partly arisen from our lack of infor- 
mation about the redshift distribution of very faint field galaxies. Most of the quantitative analyses 
undertaken adopt pure-luminosity evolution models to describe N{z) for the faint galaxy samples. 
These models, even in their milder forms, are in conflict with the observed redshift distribution in 
faint field surveys (e.g. Glazebrook et al. 1994) and their applicability at even fainter magnitudes 
must, therefore, be in doubt. At the depths probed by the current samples use to estimate u>(6), it 
is likely that even 10-m class telescopes will be unable to provide complete spectroscopic samples. 
Nevertheless, some distance information is available for these very faint sources. 

Firstly, spectroscopic redshifts have been measured for a small, but growing, sample of giant 
gravitationally lensed arcs. These arcs are highly distorted images of serendipitously placed galax- 
ies seen through the cores of rich, moderate redshift clusters. They should, therefore, represent an 
unbiased sample of high-redshift faint field galaxies. This has been confirmed by Smail et al. (1993), 
who determined optical and optical-infrared colours of arcs and found that the colours were repre- 
sentative of the bulk of the faint field population. The current sample of giant arcs has a modest 
median redshift of {z) ~ 1 with intrinsic source magnitudes of B ~ 25-26. 

Statistical information is available on the distances of samples much larger than that of the 
giant arcs from analyses of weak lensing by rich clusters of galaxies (Smail et al. 1994; Kneib 
et al. 1994). Both of these studies indicate that the bulk of the B ~ 26-27 population lies at 
redshifts z ~ 1. In particular, the analysis of an / selected sample by Smail et al. (1994) gives a 
preferred redshift distribution peaked around (z) ~ 0.8, similar to the distribution of the giant arc 
sample. The general shape of the redshift distributions indicated by the various lensing analyses 
are close to their respective NE model predictions. While the NE model is physically implausible 
at the distances and, hence, look-back times probed at these faint limits, the form of the predicted 
distribution is similar to that expected from the Burst model of Broadhurst, Ellis & Shanks (1988) 
(see also Babul &: Rees 1992). This model effectively steepens the faint end slope at moderate 
redshifts by invoking short periods of intense star formation in intrinsically faint galaxies^] to boost 
the faint number counts, while retaining the general form of the NE redshift distribution (see also 
McGaugh 1994). The faint counts are thus dominated by a population of blue, star-forming dwarf 
galaxies at modest redshifts. We therefore adopt the form of the NE redshift distribution as the 
N(z) for our standard model in the analysis below. 

Another currently fashionable galaxy evolution model, which predicts an N(z) for the faint field 
population close to that observed, is the Merger model of Broadhurst, Ellis & Glazebrook (1992). 
Unfortunately, the clustering evolution of the merging galaxy population in this model is likely to 
be complex and little theoretical work has been done on predicting this evolution. We will therefore 
not compare our observations with this model in the discussion below. 

RSMF claim that at their very faintest limits (B ~ 25) oj{6) reaches a lower limit and then 

In the discussion that follows we will use the term 'dwarf for these intrinsically low-luminosity and low surface 
brightness (LSB) galaxies; this is used with no implication for the morphological nature of the systems. 
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flattens out. Such a flattening, if it is confirmed, might indicate the presence of a high-redshift cut- 
off in the distributions of faint galaxies, either due to galaxy formation or the Lyman limit entering 
the bandpass for high-redshift sources. The existence of such a large population of high-redshift 
galaxies would not be consistent with the lensing-derived N(z). The simplest method to distinguish 
between these possibilities is to repeat this measurement in a redder passband. Therefore, we have 
undertaken a study of the evolution of the clustering strength of galaxies in a sample of faint galaxies 
selected in a red passband (Gunn r) to an effective magnitude limit deeper than all previous studies 
(r ~ 26). By combining these observations with the latest constraints on the distance to the faint 
galaxy population from lensing studies, we hope to gain new insight into the clustering evolution 
of distant galaxies. Our dataset also differs from those used previously in that it has both good 
seeing and high spatial sampling, resulting in negligible confusion at our faintest limits. 



2 OBSERVATIONS 

The observations analysed here were originally taken to select targets for a deep spectroscopic 
survey with the 10-m Keck-I telescope. The data have also been used to study the coherent 
distortion of faint galaxy images resulting from weak gravitational lensing by large-scale structure 
(Mould et al. 1994). The dataset and its reduction to a catalogue of detected objects is detailed in 
Mould et al. (1994) and here we present only a brief outline. 

The final dataset consists of a total of 24.0 ksec integration in Gunn r and 6.0 ksec in Gunn 
g on a single blank field. This was all taken under good conditions (seeing 0.7-0.9 arcsec) with 
the COSMIC imaging spectrograph on the 5-m Hale telescope. In direct imaging mode, COSMIC 
has a 9.5x9.5 arcmin field with 0.28 arcsec/pixel sampling. The final stacked r image has a la 
surface brightness limit of \i r = 28.8 mag arcsec -2 , seeing of 0.87 arcsec FWHM, and a total area 
of 90.1 sq. arcmin. The object catalogue created from this frame using the FOCAS image analysis 
package (Valdes 1982) contains ~6600 objects brighter than the 80% completeness limit of r = 26.2. 
Adopting an extremely conservative magnitude limit of r < 26.0, where the detections are roughly 
97% complete, we obtain a cumulative surface density of 71.8 galaxies per sq. arcmin (5-7 times 
the NE prediction, depending upon Qq). To calculate the equivalent R limit we use the typical 
galaxy colour at this depth and the photometric conversion of Kent (1985) giving R ~ r — 0.55 
and a limit of R < 25.5. The g-band data are deep enough to provide g — r colours accurate to 
A(g — r) ~ 0.2 at this limit. Our magnitude definition follows that of Lilly et al. (1991), adopting 
isophotal magnitudes until an object's isophotal diameter shrinks below 3 arcsec, at which point 
photometry within a fixed 3 arcsec aperture is used. The catalogue of galaxies with r < 26.0 is 
used as the basic data for the analysis detailed below. 

A mask frame defining areas around bright objects and along the frame boundaries, where the 
efficiency of galaxy detection is lower than average, was also constructed. Approximately 5% of the 
total area of the frame is masked out and these regions are excluded in the analysis below. 



3 ANALYSIS 

3.1 Estimation of uj{6) 

We estimate to(0) for the galaxies using the direct pair-counting method proposed by Landy & 
Szalay (1993), 

where DD, DR, and RR are the number of distinct data-data, data-random, and random-random 
pairs (appropriately scaled by the number of data and random points) in a given angular separation 



4 



bin. To estimate uj(9) we use 50,000 random points over the geometric area covered by the data 
catalogue, taking into account the regions which are masked out due to the presence of bright 
objects or the frame edge. The mean Cj{9) obtained from 50 independent sets of random points is 
computed for angular scales, 17" < 9 < 427", valid for our dataset (0.14 to 3.5 Mpc at the median 
redshift of our adopted N(z)). Errors for Cj{9) are estimated using 50 bootstrap resamplings of the 
data (e.g. Barrow, Bhavsar, & Sonoda 1984). 

As a consistency check of uj(9) obtained using the direct estimator above, we have applied a 
counts-in-cells estimator 

arm- Ml i 

where iVj and Nj are the number of galaxies found in cells i and j, and the brackets denote an 
average over pairs of cells separated by an angle 9 ± 59. The galaxies were binned in 5" x 5" cells 
and Co{9) computed using equation (3) above. Excellent agreement was found between the results 
of the direct and counts-in-cells estimators. 

3.2 A w and the Integral Constraint 

In the calculation of oj{9) the observed number density of galaxies on the CCD frame (to a 
given magnitude limit) is used to estimate the true mean density of galaxies at that magnitude 
limit. The small area of the frame results in a bias in the estimate of uj(9) due to the well-known 
"integral constraint" which has the effect of reducing Qj{9) by an amount 

C = ~kJ I u{ V) d ^ d ^ ( 4 ) 

where the integrals are performed over the total solid angle, Q, of the regions of the field not 
excluded by the detection mask. 

Assuming a power law form for the two-point correlation function, u){9) = A L0 9~ S with 5 = 0.8, 
as suggested by previous brighter surveys, we find C = 0.0146A,, for our field geometry and galaxy 
detection mask. 

In order to compare estimates of uj(9) from catalogues with differing depths, it is common to 
quote the value of u>(9 = 1°) after adopting a power law form for uj(9) with 5 = 0.8 and correcting 
A w for both the integral constraint and the dilution effect of faint stars which have not been removed 
from the object catalogues. The corrected amplitude is then 

.corr _ ( N ohj \ ,IC (r\ 

where A Q bj is the total number of objects used in the determination of u)(9), N s is the number of 
contaminating stars (estimated from the model of Bahcall Sz Soneira 1980), and A]p is the amplitude 
of the best-fit power law with 6 = 0.8, corrected for the integral constraint. Table 1 summarizes our 
results for these quantities, where the value of A™ rr is appropriate for 9 in units of arcseconds. The 
error on A c ° rr has been estimated from Monte Carlo simulations of fields populated with galaxies 
distributed with a known u>{9) of the form A UJ 9~ ' 8 . These simulated fields were constructed using 
the iterative tree method of Soneira & Peebles (1978). The resulting fractional error in A^ rises 
from 10% for our brightest bin to 20% for the faintest. The stellar contamination correction from 
the Bahcall & Soneira (1980) model is highly uncertain; however, a factor of 2 increase the number 
of stars in our faintest bins amounts to only a 1 <7 increase in A u . 
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3.3 Expected Correlations 



From local bright, optically-selected galaxy surveys, the two-point spatial correlation function, 
£(r), is well-approximated by a power law of the form ^(r) = (r/r ) -7 , 7 ~ 1.8 and r ~ 5.5/i _1 Mpc, 
for 10fo _1 kpc <; r & 10/i _1 Mpc (/i is the present value of the Hubble parameter in units of 100 kms 
sec -1 Mpc -1 ). The angular correlation function is related to the spatial correlation function by an 
integral equation (see, for example, Peebles 1980). Parameterizing the evolution of £(r) by 



r 



-7 

-(3+6) 



Z(r,z) = [-) (l + z)-^\ (6) 



the relation between uj{9) and £(r, z) for small angles is 
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Here x is the coordinate distance at redshift z, r is the proper coordinate, dN/dz is the number of 
galaxies per unit redshift, and the metric is 

ds 2 = c 2 dt 2 - a 2 [dx 2 /F(x) 2 + x 2 d6 2 + x 2 sin 2 ## 2 ]. 

For a power law parameterisation with constant 7, the evolution of £(r) with redshift is given by 
e. Assuming 7 = 1.8, linear theory predicts e = 0.8, while clustering fixed in proper coordinates 
yields e = 0.0, and clustering fixed in comoving coordinates yields e = —1.2. 

The small-scale clustering of galaxies is a highly non-linear process, and it is not at all obvious 
what value of e to choose in order to compare model predictions of u>(9) to observations. Neither 
is it apparent that the assumption of a unique power law form for £(r, z) is appropriate. Relatively 
few numerical studies of the evolution of £(r) from redshifts of order z ~ 3-4 to the present have 
been performed. Yoshii, Peterson, & Takahara (1993) investigated the evolution of £(r) of "density 
peak tracers" in N-body simulations of cold dark matter (CDM) universes from z = 2 to z = 
and concluded that a flat CDM model is in agreement with the low clustering amplitude of faint 
galaxies observed by EBKTG. 

Brainerd & Villumsen (1994) used a large N-body simulation to investigate the evolution of £(r) 
of individual "dark matter halos" in a flat CDM universe from z = 5 to z = and found that the 
evolution of £(r) for these objects depended strongly on the parameters used to classify groups of 
particles as "halos" (eg. mass and overdensity) . The correlation functions of moderate overdensity 
halos (Sp/p ~ 250) and high-mass, high-overdensity halos {5p/p ~ 2000) were well-described by 
a power law whose index, 7, remained constant over the course of the simulation. However, 7 
for low-mass, high-overdensity halos increased continuously. In all cases, the "first generation" 
of halos was highly clustered, in agreement with the predictions of biased galaxy formation (eg. 
Bardeen et al. 1986). Over the epoch of halo formation, however, the distribution of moderate 
overdensity halos became less clustered than that of the first generation, while the very overdense 
halos remained at the same clustering strength. For those halos whose £(r) was well described 
by a unique 7, two general patterns in the evolution of £(r) were observed after all halos had 
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formed: e = 0.0 ± 0.2 (low-mass, moderate-overdensity halos) and e = — 1.2 (high-mass, moderate 
overdensity halos and high-mass, high-overdensity halos). In the following analysis we adopt these 
values of e, along with the linear theory prediction e = 0.8, as representative of likely evolution in 
£(r). We do note, however, that more rapid rates of clustering evolution have been observed in 
N-body simulations (eg. Melott 1992, Davis et al. 1985). 

4 RESULTS AND DISCUSSION 

Figure 1 shows Cj{6) for our three faintest samples together with the best-fit power laws, A^O' 6 , 
all of which have indices of <5 ~ 0.8. This indicates that the #~°- 8 power law form for Co{6) measured 
at bright magnitudes continues all the way to our faintest sample. As observed by previous inves- 
tigators, the amplitude of to(0) decreases with the effective depth of the catalogue. We illustrate 
this with Figure 2, in which we have plotted d>(9) extrapolated to 6 = 1° (using 6j(8) oc 8~ - 8 ) as 
a function of limiting R magnitude for a number of red-selected samples. It is readily apparent 
that our results extend the trend of weaker correlation at fainter depths, with no indication of any 
flattening of this decline to our faintest limits. A weighted linear least squares fit to the data in 
Figure 2 yields u(9 = 1°) oc ^-0- 2 7±o.oi_ yyhile these results are from only a single field, they are 
in qualitative agreement with preliminary results from a larger study by Bernstein & collaborators 
(Bernstein, priv. comm.). 

The distribution of object colours on our frame radically changes at r ~ 23, when the median 
colour shifts from g — r ~ 1 to (/- r~ 0.1-0.2. Using these g — r colours we, therefore, split the 
sample by colour at g — r = 0.3 and determine the relative clustering strengths of the "red" and 
"blue" galaxies. To within the measurement errors, there is no apparent difference in u{9) between 
the two subsamples. 

We compare our observed amplitudes of lu(8) with that predicted by our standard model, a 
single population of galaxies which have a present-day correlation length typical of local bright 
galaxy samples. For simplicity this population is assumed to have the same mix of galaxy types 
visible at each epoch. To make this comparison we calculate the appropriate NE N(z) for a surface 
brightness selected sample of galaxies in each of our magnitude intervals (King & Ellis 1985) and, 
for convenience, we parameterise the shape of these predicted redshift distributions by (EBKTG) 



(It should be noted that the predicted correlation amplitude is dependent upon only the shape of 
N(z) and not its normalisation.) Adopting 7 = 1.8, ro = 5.5/i -1 Mpc, and choosing a value of e, 
we calculate for different magnitude ranges using equation (7) above. The model parameters zq 
and (3 are summarized in Table 2 for the cases of (1) a flat universe with £Iq = 1 and (2) an open 
universe with Qq = 0.2. 

Figure 3 shows our observed variation in A" m with limiting magnitude and model predictions 
using N{z) from Table 2 with the choices of e = —1.2, 0.0, 0.8 and fio = 1-0, 0.2. These models 
predict an amplitude that is about an order of magnitude greater than that observed. Clearly at 
least one of the assumptions in the models is incorrect. We first investigate the individual variation 
in each of the main model parameters (e, N(z), and ro) necessary to match our observations. 

If we allow only e to vary in the models, then we find it must be of order 6 or 7 for the 
predicted amplitude to match that observed (e s tan in Table 2), implying a very rapid evolution in 
the correlation function: £(z) ~ (1 + z)~ 7 ■ Such extreme evolution in the clustering of the entire 
population would have been easily detected by Cole et al. (1994) and, so, we prefer to search for 
another explanation. 




(8) 
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One possibility is to relax our limits on the redshift distributions. By distributing the galaxies 
over both a wider range of redshifts and out to higher redshifts, it is possible to weaken the predicted 
clustering amplitude substantially. Using the measured correlation amplitude and assuming e = 0.8, 
we determine the minimum values of zq and (3 that are consistent with the 95% confidence level 
upper limit on the measured amplitude. This choice of e allows for clustering that evolves at the rate 
predicted by linear theory and, amongst our 3 standard choices for e, yields the lowest amplitude. 
To further constrain the possible distributions, in particular the value of /3, we require that N(z) 
comply with the limits on the fraction of galaxies with z > 3 from the study of Guhathakurta, 
Tyson & Majewski (1990). The minimum allowable parameters of the model N(z) are given in 
Table 3 and indicate that by putting the galaxies at a higher median redshift, the growth of £(r) 
required to match our observed amplitude of Co{9) is much less extreme than that predicted using the 
N{z) suggested by the lensing studies. However, to achieve the necessary dilution of the clustering 
strength we require the median redshifts in our sample to be roughly 3-4 times the values obtained 
from the lensing analyses. This seems unlikely, although such a distribution might be compatible 
if the scale sizes of faint galaxies decrease very strongly with redshift. 

In our standard model we made the simplifying assumption that at each epoch we would observe 
the same mix of morphological types. This is not valid in our passband at the depths probed, even 
for the NE model. The 4000 A break moves through the r-band at a redshift of z ~ 0.7, and 
this is the median redshift proposed for our faintest samples. The presence of strong 4000 A 
breaks in early-type galaxies compared to late-types, in the absence of luminosity evolution, causes 
our deepest samples to be dominated by late-type spiral and irregular galaxies. At r ~ 20 the 
mix of types in the NE model is (E/SO/Sab, Sbc/Scd/Sdm): (98%, 2%). However, by r ~ 26 
this has changed to (37%, 63%). In the Burst model this shift is even more extreme, as the 
'excess' population in the faint counts consists exclusively of bursting gas-rich late-type spirals. This 
preeminence of moderate-redshift late-type systems in our catalogue is supported by the observed 
colour distribution, which peaks at g — r ~ 0.1-0.2, typical of late-type spirals and irregulars at 
z ~ 1. If these late- type galaxies are significantly less clustered than our adopted local sample, we 
would observe an amplitude for the correlation function which gradually falls below the predictions 
of our standard model at fainter magnitudes. 

There exist local late-type galaxy populations that exhibit weaker clustering than the bright 
optical galaxy samples used in this comparison (e.g. dwarf galaxies, LSB galaxies, IRAS galaxies). 
These could serve as a natural mechanism for diluting the clustering strength of the faint galaxy 
samples through a progressive domination by weakly clustered, late- type dwarf galaxies. We can 
predict the necessary correlation length of these galaxies must have to agree with our observations. 
Summarised in Table 4 are the correlation lengths, ro, required of this population in order to 
match the 95% confidence upper limits on our observed amplitudes. The choice of an open or 
closed universe makes little difference in the required correlation length, and for all three values of 
e, ro is less than half the present-day correlation length of bright, optically-selected galaxies. 

By a similar argument Bernstein et al. (1994) were able to construct a viable model for their 
observations by adopting the correlation function of IRAS galaxies for the dominant blue population 
in their B ~ 22 sample. They then required only modest clustering evolution (e ~ 0) for the bulk 
of their z ~ 0.3 population. Shown in Table 2 is the clustering evolution needed to reproduce our 
values of A u using a population described by the IRAS correlation function of SRRL (7 ~ —1.6 and 
ro ~ 3.8/i _1 Mpc). Under these assumptions our correlation function has to evolve as eiRAS ~ 3, 
which is slower than our standard model predicts, but considerably faster than the theoretical 
predictions. 

We would claim that the IRAS correlation function is probably not the correct choice to de- 
scribe the faint galaxy population. A more representative local sample would be either dwarf or 
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LSB galaxies (see McGaugh 1994), especially within the framework of the Burst model described 
above. Two recent, large studies have measured correlation functions for similar classes of galaxies. 
Santiago & da Costa (1990) used volume-limited samples from the new Southern Sky Redshift 
Survey (SSRS) to study the clustering of LSBs, while Thuan et al. (1991) analysed a mixed sample 
of dwarf and LSB galaxies from the UGC. The more interesting of these is the SSRS study which 
finds a correlation length for LSB galaxies of vq ~ 2.3-2.7/i _1 Mpc, depending upon the exact 
sample definition. These local values are somewhat uncertain due to the small samples available, 
however, the range of values is in reasonable agreement with those quoted in Table 4 for clustering 
growth of e ~ 0.8. It therefore appears that we can achieve a sufficiently low amplitude for lo(9) at 
faint limits by postulating a gradual transition to LSB, late-type dominated samples via differential 
if -corrections, coupled with the intrinsically weaker clustering of these galaxies. 

5 CONCLUSIONS 

We have measured the angular correlation function of galaxies to an apparent magnitude limit of 
r < 26. The correlation amplitude, A^, for our faintest sample is extremely small, but significantly 
non-zero. This detection is approximately an order of magnitude weaker than the predictions of 
a theoretically and observationally motivated standard model. The amplitude of the correlation 
function decreases continuously to our faintest limits, showing no evidence for any flattening (c.f. 
RSMF). We find no difference between the clustering strength of the blue and red galaxies within 
our sample. 

There are three main parameters in our standard model which can be varied in order to obtain 
agreement with the observations: the evolution in £(r) (parameterised by e), the redshift distri- 
bution of the sample, N(z), and the correlation length, rQ. Forcing the observed low correlation 
amplitude to result from a variation of either e or N(z) results in two unpalatable solutions: ex- 
tremely rapid (and unexpected) evolution in £(r) or a redshift distribution in which the galaxies 
are distributed at median redshifts 3-4 times that expected from gravitational lensing studies of 
rich clusters. 

We prefer a more natural explanation for the low clustering strength observed, resulting from 
the tendency of deep field samples to be dominated by moderate-redshift, blue star-forming spiral 
and irregular galaxies. This arises from the different if -corrections for early- and late-type systems. 
However, this effect is magnified by the overwhelming prevalence of blue late-type systems in faint 
field counts, as illustrated by the trend to bluer colours at fainter magnitudes. In the Burst model 
of Broadhurst, Ellis & Shanks (1990) these late-type systems are associated with star bursting LSB 
and dwarf galaxies. Assuming modest clustering growth we derive a correlation length for these late- 
type systems of ro ~ 2h~ l Mpc. This is significantly shorter than the equivalent value for either 
local bright optically-selected galaxies, or IRAS-selected galaxies. However, it is in reasonable 
agreement with the correlation lengths measured for LSB galaxies from the recent SRSS (Santiago 

6 da Costa 1990). 

We find reasonable support for a simple model of faint galaxy evolution where luminous galaxies 
are distributed out to modest redshifts (z ~ 1) and the high galaxy surface density is made up of 
a dominant co-eval population of bursting LSB galaxies. The clustering of this dominant, bursting 
population can be explained by plausible theoretical clustering evolution of the locally observed 
samples. 

A simple observational test of our model for the clustering of faint galaxies would be a correlation 
analysis of a deep infra-red (K) selected sample, which should show significantly stronger clustering 
than we observe. This results from the dominance, out to high redshifts, of the same early-type 
galaxies which define the local bright optically-selected samples. 
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FIGURES 

Figure 1: The two-point angular correlation function of galaxies in 3 magnitude intervals: (a) 
20.0 < r < 25.0, (b) 20.0 < r < 25.5, (c) 20.0 < r < 26.0). No correction for the integral constraint 
has been made. Error bars are estimated from 50 bootstrap resamplings of the data. The best-fit 
power laws of the form uj(9) = A6~ - 8 are indicated by the dotted lines. Note that different vertical 
scales are used in each panel. 

Figure 2: Corrected ui{9) extrapolated 1 degree as function of limiting R magnitude for a number 
of red selected samples. The magnitude limits for the CJB study have been converted to R using 
the transformations given by Yoshii, Peterson & Takahara (1993) and RSMF. This figure is adapted 
from Roche et al. (1993). 

Figure 3: Observed correlation function amplitude, ^4™ rr , as a function of limiting r magnitude 
(squares). Error bars indicate 95% confidence limits. The predicted amplitude as a function of 
limiting magnitude, e, and f^o (assuming ro = 5.5/i _1 Mpc and 7 = 1.8) is indicated by lines (solid 
for Qq = 1.0; dashed for Qq = 0.2). Also shown (dotted line) is the predicted A C ° TT assuming 
r = 2.0/i- 1 Mpc, Q = 1.0, 7 = 1.8, and e = 0.8. 
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Table 1: Observed correlation function parameters 



mag. 


interval 


A^obs 


( Wobs y 

\N oha -N s ) 


4corr 


u(9 = i°) 




20.0 < 


r < 24.5 


2170 


1.439 


0.37±0.04 


(5.3 ± 0.6) x 10" 


-4 


20.0 < 


r < 25.0 


3103 


1.329 


0.31±0.04 


(4.4 ±0.6) x 10" 


-4 


20.0 < 


r < 25.5 


4292 


1.258 


0.20±0.03 


(2.9 ±0.4) x 10" 


-4 


20.0 < 


r < 26.0 


5730 


1.206 


0.11±0.02 


(1.6 ±0.3) x 10" 


-4 



Table 2: Parameters for N(z) assuming a NE model and corresponding evolution of to(9) 









= 1.0 








= 0.2 




mag. interval 


zo 


P 


e stan 


eiRAS 


zo 


P 


e stan 


eiRAS 


20.0 < r < 24.5 


0.49 


2.5 


7 


3 


0.45 


2.4 


7 


3 


20.0 < r < 25.0 


0.58 


2.7 


6 


3 


0.48 


2.4 


7 


3 


20.0 < r < 25.5 


0.65 


2.7 


6 


3 


0.55 


2.5 


7 


3 


20.0 < r < 26.0 


0.75 


2.8 


7 


4 


0.64 


2.7 


7 


4 



Table 3: Shallowest N(z) consistent with the 95% confidence level upper limit on lj(0) 





^0 = 


1.0 


^0 = 


0.2 


mag. interval 


zo 


P 


zo 





20.0 < r < 24.5 


1.4 


3 


1.4 


6 


20.0 < r < 25.0 


1.8 


4 


1.6 


6 


20.0 < r < 25.5 


2.4 


6 


1.8 


6 


20.0 < r < 26.0 






2.5 


6 



Table 4: Correlation length, vq (h 1 Mpc), required to match A c °™ assuming a NE model N(z) 
and S = 0.8 (0.6) 







Q = 1.0 






fi = 0.2 




mag. interval 


e = -1.2 


e = 0.0 


e = 0.8 


e = -1.2 


e = 0.0 


e = 0.8 


20.0 < r < 24.5 
20.0 < r < 25.0 
20.0 < r < 25.5 
20.0 < r < 26.0 


1.4 (1.7) 
1.4 (1.7) 
1.2 (1.4) 
0.9 (1.0) 


1.8 (2.2) 
1.8 (2.3) 
1.6 (2.0) 
1.3 (1.5) 


2.1 (2.6) 

2.2 (2.8) 
2.0 (2.4) 
1.6 (1.9) 


1.5 (1.8) 
1.5 (1.8) 
1.3 (1.5) 
1.0 (1.1) 


1.9 (2.3) 
1.8 (2.3) 
1.7 (2.0) 
1.3 (1.5) 


2.2 (2.7) 
2.2 (2.8) 
2.0 (2.4) 
1.6 (1.9) 
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Fig. 1 



.03 



"i i i i i i i i i i i i i i i i i i i i i i i i i r 



.02 — 



.01 



(a) 



I 



J I I I I I I I I I I I I I I I I I I I I I I I I L 

1.4 1.6 1.8 2 2.2 2.4 



.015 



3 




1.4 1.6 1.8 2 2.2 2.4 

"i i I i i i I i i i I i i i I i i i I i i i I i i r 



.01 



.005 



(c) 



I 



J I I I I I I I I I I I I I I I I I I I I I I I I L 

1.4 1.6 1.8 2 2.2 2.4 



iog 10 (e") 



Fig. 2 
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